clc;
clear;
load x_actual_pcc.txt;
load pcc_esti_xekf.txt;
load pcc_esti_xekf1.txt;
load pcc_esti_xekf2.txt;
% load pcc_esti_xekf3.txt;

actual = x_actual_pcc(:,[68,69,73,80,85,88]);
esti = pcc_esti_xekf(:,[68,69,73,80,85,88]); % optimal
esti_1 = pcc_esti_xekf1(:,[68,69,73,80,85,88]);
esti_2 = pcc_esti_xekf2(:,[68,69,73,80,85,88]);
% esti_3 = pcc_esti_xekf3(:,69:88);

t = 1:2:180;
figure(1)
for i = 1:6
    subplot(2,3,i)
    h1 = plot(t,actual(:,i),'k','Linewidth',1); hold on; 
    h2 = plot(t,esti(:,i),'color',[0/255 161/255 59/255],'Linewidth',1); hold on;
    h3 = plot(t,esti_1(:,i),'m','Linewidth',1); hold on;
    h4 = plot(t,esti_2(:,i),'r','Linewidth',1); hold on;
%     h5 = plot(t,esti_3(:,i),'b','Linewidth',1); hold on;
    legend([h1,h2,h3,h4],'Actual','Case 1','Case 2','Case 3')
    xlabel('t/min');
    xlim( [ 0, 182 ] );
end

subplot(2,3,1)
ylabel('$C_3^G(CO_2)$',  'interpreter','latex');
subplot(2,3,2)
ylabel('$C_4^G(CO_2)$',  'interpreter','latex');
subplot(2,3,3)
ylabel('$C_3^G(MEA)$',  'interpreter','latex');
subplot(2,3,4)
ylabel('$C_5^G(H_2O)$',  'interpreter','latex');
subplot(2,3,5)
ylabel('$T_5^G$',  'interpreter','latex');
subplot(2,3,6)
ylabel('$T_{reb}^G$',  'interpreter','latex');


% ylabel('CPI');
